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Abstract. Chiral effective field theory (EFT) provides a systematic and controlled approach 
to low-energy nuclear physics. Here, we use chiral EFT to calculate low-energy weak Gamow- 
Teller transitions. We put special emphasis on the role of two-body (2b) weak currents within 
the nucleus, and discuss their applications in predicting physical observables. 



1. Introduction 

Electro-weak nuclear reactions at low energies, which are the subject of this contribution, have 
an important role in many physical scenarios. They are the microscopic motor of astrophysical 
phenomena, such as solar fusion, core collapse supernovae and the nucleosynthesis within. In 
addition, these reactions are used to study nuclear structure and dynamics. Moreover, low- 
energy /3-decays and weak transitions provide a window into high-energy physics. For example, 
superallowed decays allow high precision tests of the Standard Model, and 0^/3/3 decays probe the 
nature of neutrinos, their hierarchy and mass. High accuracy theoretical, as well as experimental 
understanding of the nuclear response is a key ingredient needed to accomplish these tasks. 

However, the non-perturbative character of QCD at low energies makes a theoretical 
calculation of the nuclear response extremely challenging. In the last two decades, chiral 
effective field theory of QCD at low energies was developed, enabling a construction of a nuclear 
Lagrangian that is consistent with QCD. Chiral EFT provides a systematic basis for nuclear 
forces and consistent electroweak currents [TJ [2] , where pion couplings contribute both to the 
electroweak axial current and to nuclear interactions. This is already seen at leading order: 
the axial constant, gA, determines the axial one-body (lb) current and the one-pion-exchange 
potential. Two-body (2b) currents, also known as meson-exchange currents, enter at higher 
order, just like three-nucleon (3N) forces pp. As shown in Fig. [TJ the leading axial contributions 
are due to long-range one-pion-exchange and short-range parts [2], with couplings C3, C4 and cd, 
which also enter the leading 3N (and subleading NN) forces [U [3] . 

In this contribution we present some recent studies of chiral currents and their consistency 
with nuclear forces. We emphasize the importance of using chiral 2b currents in weak Gamow- 
Teller transitions in nuclei. 

2. Gamow- Teller transitions 

At low energies, the coupling of weak probes to a nucleus is given by the current-current 
interaction, H\y = ^| J d 3 r e~ %v ' v ini + h.c, where Gp is the Fermi constant, p the 
momentum transferred from nucleons to leptons, and jh^ the leptonic current of an electron 
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Figure 1. In chiral EFT, Feynman diagrams of chiral 2b currents (left) correspond to 3N force 
contributions (right) that depend on the same low energy constants. 

coupled to a left-handed electron neutrino. This approximation creates a separation between 
the leptons and baryons, which is a great simplification of the problem. It suggests that only the 
symmetry of the interaction is needed to determine the structure of the current. For example, 
any axial probe, regardless of its fundamental origin, would couple to the same current within the 
nucleus, which is the Nother current derived from the axial symmetry of the chiral Lagrangian, 
which we will consider up to next-to-next-to-next-to leading order (N 3 LO) [2]. 

The nuclear current jft is organized in an expansion in powers of momentum Q ~ m n over 
a breakdown scale Ab ~ 500 MeV. To order Q 2 (and also Q s in this counting), the lb current, 
= Eii T r[^° J °ib - 6fikj i,ib\ 6 ( r - r *)> has temporal and spatial parts in momentum 
space [2]: 

t0 / 2\ ( 2\ P-<T« , 2 E(p-(Ti) 

Ji,ib(P) = 9v{p)-9a^ — +gp(p) — 7, , I 1 ) 
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JiMP) = 9a(p )<7i- gp(p ) — ^ — + i(9m + gv) - gv ^ , (2) 

where E = E t - E[, p = Pi - p[, and P = Pi + p[. g v (p 2 ), gA(p 2 ), gp(p 2 ), and gM(p 2 ) are the 
vector (V), axial (A), pseudo-scalar (P), and magnetic (M) couplings. In chiral EFT, the p 
dependence is due to loop corrections and pion propagators, to order Q 2 : gv,A(p 2 ) = gv,A (1 — 
2-g-), with g v = l,A v = 850 MeV, A A = 2^/r A = 1040 MeV; g P {p 2 ) = ^p^-4^(p 2 )^ 

and qm = fJ'p — P-n = 3.70, with pion decay constant F n = 92.4 MeV, m n = 138.04 MeV, and 
g npn = 13.05 [I]. 

At leading order Q°, only the momentum- independent g A and gy terms contribute. They 
give rise to low- momentum-transfer Gamow- Teller and Fermi decays. The former has the known 
form,^| L o = ^Ef=i^r. 

Among the Q 2 terms, form-factor-type (ff ) contributions and the gp part of J^ib dominate. 
The remaining Q 2 terms are odd under parity, so they require either a P-wave lepton or another 
odd-parity term to connect to + states, common for ground states. In such cases, the P and 
E terms in Eqs. ([T]) and ^ can be neglected, and only the term with the large gu + gv = 4.70 
is kept. 

At order Q 3 , 2b currents enter in chiral EFT [2]. These include vector spatial, axial temporal, 
and axial spatial parta^ The first two are odd under parity, and therefore can be neglected 
in couplings to + ground states. In addition, within chiral EFT the vector spatial current 
is conserved, satisfying the CVC hypothesis, thus at low-energies their contribution can be 
calculated via the Siegert theorem. As a result, the dominant weak 2b currents include an axial 
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Table 1. Calculated 3 H, 3 He and 4 He g.s. energies (MeV) and point-proton root-mean-squared 
radii (fm) 
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Equation ^ includes contributions from the one-pion-exchange C3 , C4 parts and from the short- 
range couplings d\,d,2, where due to the Pauli principle only the combination d% + 2d 2 = 
c d/{9A^-x) en ters [with A x = 700 MeV], see Fig [TJ This result is extremely interesting as 
it connects the low-energy constants (LECs) determining the strength of the meson exchange 
current (MEC) contact term to cjj, which is one of the two unknown low-energy constants, cd 
and ce, that determine the 3N force up to N 3 LO, and do not contribute to 2N force. This 
relation shows that one can use a weak observable, e.g., an empirical extraction of (E^ 1 
from a decay rate, as a constraint for the determination of cd and ce- 
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3. Using triton /3-decay to determine the 3N force 

We use the relation between the short range contribution to the current and the 3N force to 
fully constrain the Lagrangian in systems with up to 3-body, and predict 4 He properties. The 
calculation steps are as follows: (i) calculate the 3 H and 3 He g.s. wave functions by solving 
the Schrodinger equation for three nucleons interacting via the chiral iViV potential at N 3 LO 
of Ref. [5] and the NNN interaction at N 2 LO [6]; (ii) determine for which cd values along the 
trajectory the calculated reduced matrix element of the E± operator (at N 3 LO) reproduces 
(E^) emp extracted from 3 H /3-decay. We find that these two constraints are uncorrelated, 
producing a stringent determination of the LECs: cd = —0.2 ±0.1 and ce = —0.205 ±0.015 |3j. 

With this calibration of cd and ce, for this potential, in principle, any other calculation is 
a prediction! In Table [T] a collection of A = 3 and 4 data is given, obtained with and without 
inclusion of the NNN force for cd = -0.2 (c e = -0.205). The NN+NNN predictions for the 
4 He are in good agreement with measurement. 

To this important conclusion, one can add the results of few additional tests, appearing in 
Fig. [2j that are aimed to probe the axial correlations in the nucleus, by analyzing the sensitivity 
of the triton half life to NNN force and/or MEC: 

• The fundamental importance of the axial two-body currents. By suppressing the MEC, in 
the whole investigated cd-ce range, the calculations underpredict {Ef) emp by about 2%. 
Within the MEC, one finds a large cancellation between the long-range one-pion-exchange 
term (corresponding to the left diagram in Fig. [T]), and the contact term (second diagram 
from the left Fig. [T]) . 
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Figure 2. The ratio (Ef) th /(Ef) emp using the N 3 LO NN potential [5] with and without 
the N 2 LO NNN interaction [BJ, and the axial current with and without MEC. The shaded 
area is twice the experimental uncertainty. Also shown: results for the phenomenological AV18 
potential (with A = 500 MeV imposed in the current), and for the N 3 LO AW potential of 
Epelbaum et al. [7j (with A = 450, 600 MeV, and a 700 MeV spectral-function cutoff in the 
two-pion exchange), which should be compared with the full calculation (black solid line). The 
line denoted by d,R = includes only the long range contribution to the MEC, neglecting the 
contact term contribution. Figure taken from [3]. 



• A negligible effect of the NNN force on this observable. 

• Weak dependence on the specific character of the force. 

These features, which might be unique to s-shell nuclei, entail that the determination of cd 
and ce obtained in this work is robust. 

4. Extracting the weak structure of the nucleon from fi~ capture on 3 He [8j 

In the previous section it was shown that for a weak transition in the 3-nucleon sector, the 
specific character of the force had a minor effect on the reaction rate. We extrapolate this 
understanding to the muon capture on 3 He which results in a triton, /x~ + 3 He — > + 3 H. 
We study this reaction in a hybrid approach: the calculation uses the chiral EFT based weak 
current, combined with the phenomenological nucleon-nucleon potential Argonne v\% (AV18) [9J 
augmented by the Urbana IX (UIX) [10] three nucleon force. This allows a study of the cutoff 
variation of the observable. 

[i capture on 3 He has an accurately measured rate T(^ _ + 3 He — > i/„ + 3 H) exp = 1496(4) Hz, 
i.e. a ±0.3% precision [11]. Combining with the sizable momentum transfer in ordinary muon 
capture (OMC), that enhances the effect of the induced-pseudoscalar form factor, allows its 
study. In addition, one can use OMC to study second class currents [12"] . 

The calculation process demands fixing the single contact unknown LEC, which determines 
the weak axial current. We do this by reproducing the experimental triton half-life, for various 
cutoff values for the chiral current. The results for the nuclear matrix element of the muon 
capture process show a 9% effect due to the MEC contribution, and a small 0.3% effect due 
to the cutoff dependence. It is worthwhile noting that the relative contribution of the MEC to 
this process is almost three times bigger than the MEC contribution to the triton half life. An 



extremely weak cutoff dependence shows that the essential physics is captured in the chiral EFT 
operators. 

The final prediction for the capture rate is 



where the first error is due to the Heavy Baryon chiral perturbation theory (HB^PT) cutoff, the 
second is due to uncertainties in the extrapolation of the form factors to finite momentum 
transfer, and in the choice of the specific nuclear model, the third error is related to the 
uncertainty in the triton half life, and the last error is due to theoretical uncertainty in the 
electroweak radiative corrections calculated for nuclei [13] • This sums to a total error estimate 
of about 1%, which overlaps with the uncertainty 

The conservative error estimation still allows rather interesting conclusions. First, one notices 
that the calculated capture rate agrees with the experimental measurement + 3 He — > 

+ 3 H)™^. = 1496(4) Hz. Thus, one concludes that EFT accurately predicts the capture rate. 

However, the most interesting result concerns the weak form factors of the nucleon. In order to 
constrain the induced pseudoscalar and second class form factors, we take the following approach. 
In each case, we set all the other form factors to their nominal value, and change this form factor 
in a way which keeps an overlap between the experimental rate and the theoretically allowed 
rate. The nominal value of the form factor is set to reproduce the experimental measurement. 

The resulting constraint on the induced pseudoscalar form factor is: 



in very good agreement with the HBxPT prediction of Ref. [14J. 

A second conclusion concerns the contribution of second class currents. While the current 
constraint on the axial G- parity breaking term is rather weak, this calculation puts the tightest 
limit on CVC, constraining it to ra e FsjFy = (0.5 ± 2) x 10 -4 , which is consistent with CVC. 
One notes that this is already close to the regime of CVC breaking, predicted by xPT |14j . 

A recent calculation in the framework of the previous section, i.e., taking currents and forces 
from the same chiral EFT Lagrangian, has reached similar results |15j . 

5. The in-medium behavior of the axial constant 

Surprisingly, key aspects of (3 decays remain a puzzle. In particular, when calculations of 
Gamow- Teller (GT) transitions of the spin-isospin-lowering operator g^err - are confronted with 
experiment, some degree of renormalization, or "quenching" q, of the axial coupling g A s = qg^ 
is needed. Compared to the single-nucleon value qa = 1.2695(29), the GT term seems to be 
weaker in nuclei. This was first conjectured in studies of /3-decay rates, with a typical q ~ 0.75 in 
many-body calculations |16[ I17j. In view of the significant effect on weak reaction rates, it is no 
surprise that this suppression has been the target of many theoretical works. It is also a major 
uncertainty for Ovf3f3 decay nuclear matrix elements (NMEs), which probe larger momentum 
transfers of the order of the pion mass, p ~ m n , where the renormalization could be different. 
Here we revisit this puzzle based on chiral EFT currents. 

The lightest nucleus that undergoes a /3-decay is the triton. However, the theory cannot be 
checked in the triton since its half-life is used to calibrate the strength of the MEC. The lightest 
nucleus that can provide a test to the theory is thus 6 He (J 71 " = + ), an unstable nucleus, which 
undergoes a (3 decay with a half-life t\/2 = 806.7 ± 1.5 msec to the ground state of 6 Li (J 71 " = 1 + ) 
|18j . As this is a 6-body problem, we use a soft nuclear potential, JISP16. This potential leads 
to an underbinding of about 80keV for the 3 He and 120keV for the triton, where we calibrate 
the one unknown LEC determining the current. One calculates binding energies of 28.70MeV 
for 6 He and 31.46 MeV for 6 Li. These results underestimate the experimental energies by about 



r = 1499(2)a(3) N m(5M6) rc Hz, 
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gp (q 2 = -0.954m^) = 8.13 ± 0.6 
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Figure 3. Top panel: qa plus 2b-current contributions for p = GT transitions normalized to 
gA as a function of density p. The boundaries of the box are given by q = 0.74/0.96 and 
p = 0.10, 0.12 fm~ 3 . The different curves correspond to possible calibrations of cd (see 
Ref. [IH] for details), with shaded regions for the density range. Intermediate quenching values 
would lie between these regions. The inset shows the quenching predicted by the long-range 
parts of 2b currents only (cd = 0). Bottom panel: Same as top, but as a function of momentum 
transfer p for empirical/smaller quenching q = 0.74/0.96. Figure taken from Ref. [191. 



0.5 MeV, though the difference AE = 2.76MeV differs by merely 34keV from the experimental 
value. 

Incorporating the xPT based contributions to the weak-current one can calculate the full 
6 He- 6 Li GT matrix-element, as a function of the cutoff imposed in the weak current: 

|GT( 6 He)| theo = 2.198(l) A (2) N (4) t (5) gA = 2.198 ± 0.007. (6) 

The first error is the cutoff variation dependence, the second is numerical, the third is due to 
uncertainties in the triton half- life, and the last is due to uncertainties in gA- This should be 
compared to the experimental matrix-element |GT(^He)| CX pt = 2.161 ± 0.005. Thus, the theory 
overpredicts GT by about 1.7%, which is reasonable considering the approximation in using a 
phenomenological, pure 2-body, potential. 

For heavier nuclei a microscopic calculation is still out of reach. We thus study the impact 
of chiral 2b currents in nuclei at the normal-ordered lb level by summing the second nucleon 
over occupied states in a spin and isospin symmetric reference state or core. Taking a Fermi gas 
approximation for the core and neglecting tensor-like terms (k • crk — ^k 2 <r)r~, we obtain the 
normal-ordered lb current 1191: 



eff — P 

J|,2b = -9A(TiTi 



c D 2 p 2 r , n ./l. o 1 
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where p = 2kp/(3ir 2 ) is the density of the reference state, k-p the corresponding Fermi 
momentum, and I(p, P) is due to the summation in the exchange term. The effective lb current 
Jf 2b only contributes to the GT operator and can be included as a correction to the gA{p 2 ) 
part of the lb current, Eq. ([2]). This demonstrates that chiral 2b currents naturally contribute 
to the quenching of GT transitions. Compared to light nuclei, their contributions are amplified 



because of the larger nucleon momenta. A quenching of p ~ GT transitions is predicted, as 
well as the momentum transfer dependence, as seen in Fig. [3} 

This formalism can be used to evolve the chiral EFT understanding various phenomena, 
including the prediction of 0vf3f3 decay rates. Indeed, such a calculation of the Qu/3/3 decay 
operator based on chiral EFT currents shows that 2b contributions to the nuclear matrix 
elements are significant, and range from —35% to 10% [19]. Another recent application, which 
will not be covered in this contribution, is to the scattering of a weakly interacting massive 
particle (WIMP) candidate on nuclei |20j. 

6. Summary 

We showed that chiral EFT enables parameter free calculations of electro-weak Gamow Teller 
transitions, that can be used in various fields: from constraining the nuclear forces, through a 
deeper understanding of in- medium nuclear effects, and up to predicting weak reactions for the 
use of experiments studying the limits of the Standard Model. 
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